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Abstract 

We report a preliminary branching fraction of (1.80±0.37 s t a t. ±0.23 sys t.) x 10~ 4 for the charmless 
exclusive semileptonic B + — > 7r°^ + i/ decay, where ^ can be either a muon or an electron. This 
result is based on data corresponding to an integrated luminosity of 81 fb _1 collected at the T(4S) 
resonance with the BABAR detector. The analysis uses BB events that are tagged by a B meson 
reconstructed in the semileptonic B~ — > D®l~v(X) decays, where X can be either a 7 or a 7r° from 
a D* decay. 
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1 INTRODUCTION 



Measurements at B Factories have significantly improved our knowledge of CP violation in the 
B° — B° system. In particular, the angle j3 of the Unitarity Triangle (see Figure ^) has been 
measured to a 5% accuracy from time-dependent CP asymmetries in b — > ccs decays. 

On the other hand, experimental determination of the other two angles and of the lengths 
of the two sides (with the third side normalized to unit length) have yet to achieve comparable 
precision. The uncertainty in the length of the side opposite to the angle f3 is dominated by the 
smallest CKM matrix element, \V u b\. Improved determination of \V u b\ therefore translates directly 
to a more stringent test of the Standard Model prediction. 




Figure 1: Representation of the Unitarity Triangle 



Charmless semileptonic decays of B mesons provide the best probe for | V u b \ ■ Measurements can 
be done either exclusively or inclusively, i.e., with or without specifying the hadronic final state. 
Since both approaches suffer from significant theoretical uncertainties, it is important to pursue 
both types of measurements and test their consistency. 

The exclusive B — > X u iv decay rates are related to \V u b\ through form factors (FF). In the 
simplest case of B — > nlis, the differential decay rate (assuming massless leptons) is given by 

dT(B° -> TTt+v) _ n dT(B+ -> tt^+i/) _ G F \V ub \ 2 , 2 3 

dq 2 dq 2 ~ 24vr3 j| P *> {L) 

where Gf is the Fermi constant, q 2 is the invariant-mass squared of the lepton- neutrino system and 
p n is the pion momentum in the B frame. The FF f+(q 2 ) is calculated with a variety of approaches. 
Major improvements achieved recently in the calculation of the FF, based on light-cone sum rules 
and unquenched lattice QCD [2J|3] calculations, should allow a competitive determination of \V u b\ 
using exclusive semileptonic decays. 

In this paper, we use the ISGW2 model 0j. Alternate calculations [U 12 EJ are considered to 
estimate model dependence of the result. Measurements of the branching fraction B{B — > nil/) 
have been reported by CLEO , Belle , and BABAR . The CLEO and BABAR measurements 
use neutrino reconstruction in untagged BB events; the Belle measurement uses semileptonic tags. 
BABAR has also reported a measurement of the total branching fraction B{B — > nlu) using fully- 
reconstructed hadronic tags |Hj. 

In this paper, we report a preliminary branching fraction measurement from a study of the 
B + -» tt £ + v decay, using event samples tagged by B — > D°£ v{X) decays. 4 A similar study of 
the B° — > Tt~l + v decay is reported in a separate paper [S]. 
4 Charge- conjugate modes are implied throughout this paper. 



8 



We look for combinations of a D° meson and a charged lepton (e~ or pT ) that are kinematically 
consistent with B — > D°£~V(X) decays. For each such candidate, we define the signal side as 
the tracks and neutral clusters that are not associated with the candidate. We search in the signal 
side for a signature of a B + — > 7r°£ + z/ decay. We take advantage of the simple kinematics of 
the B + — > 7r°^ + ^ process and extract the signal yield. We calculate the branching fraction using 
the signal efficiency predicted by a Monte Carlo (MC) simulation. We correct for the data-MC 
efficiency differences using control samples in which both B mesons decay to tagging modes. 

2 THE BABAR DETECTOR AND DATASET 

This measurement uses the e + e~ colliding-beam data collected with the BABAR detector at the 
PEP-II storage ring. The data sample analyzed contains 88 million e + e~ — > BB events, where BB 
stands for B + B~ or B°B°, which corresponds to an integrated luminosity of 81 fb _1 on the T(4S) 
resonance. In addition, a smaller sample (10fb _1 ) of off-resonance data recorded approximately 
40MeV below the resonance is used for background subtraction and validation purposes. 

We also use several samples of simulated BB events to evaluate the signal and background 
efficiencies. Charmless semileptonic decays B — ► X u £u are simulated as a mixture of exclusive 
channels (X u = it, r\, rf, p, and to) based on the ISGW2 model 4^ and decays to non-resonant 
hadronic states [TT] . 

3 ANALYSIS METHOD 

The event selection has been developed blind, that is, the analysis was first optimized using MC 
simulation to obtain the largest expected statistical significance of the signal yield; only then was 
it applied to data. 

3.1 Event Selection 

We search for candidate BB events in which one B meson decayed as Br -> D°£-V(X), B° -> 
D + £~V or B° — » D* + £~V, with D* + — ► Z?°7r + . The D mesons are reconstructed in the D — > K~ir + , 
K~Tr + Tr~- + , K~ir + ir°, and D + — > K~tt + - + channels. Charged B tags are used to select the 
B + — > ir®l + v decay mode. The D° candidates from charged B tags are kept if their reconstructed 
mass lies within 10a of the fitted mean of the D° mass distribution, where a is the experimental 
resolution. The signal region corresponds to |m^j° — m^o | < 3cr. We define sideband regions, 
which are used to subtract the combinatorial background, as 3a < [m^° — m^o cd | < lOcr. The 
neutral B° tags are used to reject B°B° events which are a background for the present analysis. 
In this case, a tighter D° mass selection, within 3a, is applied. The same mass window is used for 
D . The a values vary between 5MeV/c 2 and 13MeV/c 2 depending on the D decay channel. 

Candidate D* + — > D°tt + decays are required to have a mass difference, m D *+ — m D o, of ± 
5 MeV/c 2 around its expected value (145.4 MeV/c 2 ). When a D° meson is identified as being 
part of a B° — > D* + £~u(X) tag, the corresponding B~ —>■ D°£~17(X) tag is dropped. Note that 
D* — ► Dir /^ are not explicitly reconstructed. 

The D^*> candidates are combined with electrons or muons to identify a B -> D^*Uv(X) 
candidate. The lepton must have a center-of-mass momentum 5 p\ > l.OGeV/c. The charge of the 
lepton must be the same as that of the kaon from the D candidate. 

5 Variables denoted with a star (x*) are measured in the Y(4S) rest frame; others are in the laboratory frame 
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Referring to the tag-side D^l combination as the U Y" system, we calculate the cosine of #ey> 
the angle between and p y as 

_ 2E* B E Y -m%-m Y 

COS V BY — I I • (2) 

Equation © assumes that a B — > D^£u(X) decay has been correctly reconstructed, and the 
only undetected particle in the final state is the neutrino. If that is the case, cos 9 by should be 
between —1 and +1 within experimental resolution. If the tag has been incorrectly reconstructed, 
Equation ((2j) does not give a cosine of a physical angle, and cos 6 by is distributed more broadly. 
We require — 2.5 < cos 8 by < +1-1. The asymmetric range chosen for the cos 6 by cut reflects the 
fact that a large fraction of the signal events contain a B~ — > D*°£~u decay in which the soft 7r° 
or 7 was not accounted for. Such events tend to populate the cosQby distribution between —2.5 
and -1.0. 

There is often more than one remaining D&iv(X) candidate in an event. When this happens, 
we select the candidate with the smallest value of |cos#by|- The average number of D^*ku(X) 
candidates per signal event is 1.2 according to the MC simulation. If at this point the best candidate 
represents a B° decay (i.e., B° — > D^ + £~v), the event is discarded. 

From each event with a B~ — > D°£~V(X) candidate, we remove the tracks and the neutral 
clusters that compose D°£~. We then search for a B + — > 7r°£ + f candidate in the remaining part 
(signal side) of the event, which must contain an identified lepton and a pair of photons from 
7T° — > 77 that satisfy 115MeV/c 2 < m 77 < 150MeV/c 2 . The two leptons in a candidate event must 
be oppositely charged. 

In order to suppress non-Si? backgrounds, which have a more jet-like topology than BB events, 
we require the ratio R2 of the second and zeroth Fox- Wolfram moments |12| . computed using all 
charged tracks and unassociated neutral clusters in the event, to be smaller than 0.5. 

In addition, we require the momenta of the ir° and lepton to satisfy |p* 1 + |p| | > 2.6 GeV/ c. This 
cut removes more than 95% of the BB background events while retaining nearly the entire phase 
space of the B + — » k®£ + v signal. In the limit of massless leptons and 7r°s, this cut corresponds to 
the kinematic limit for the neutrino energy, E* < 2.7 GeV. 

Analogously to cos 9by, we can calculate the cosine of 9 BlT o^ the angle between p^ and p* 0£ as 

. _ 2E bKh ~ m B- ™Ih f oN 

COS#R 7r CV = • o) 

This variable, again, should be between —1 and +1 for the signal events, and distributed broadly 
for the background. We require — 1.1 < cos8 B7T og < +1.0. 

If a signal event has been correctly reconstructed in the B~ — > D°£~V and B + — > 7r°£ + ^ decays, 
no other particles should be present. In reality, such events often contain extra neutral particles. 
Some of them come from soft 7r°s and/or photons from decays of D*° or heavier charmed mesons. 
We identify the photons that may have come from D*° — > D°7r° and D°7 decays by combining 
them with the D° meson candidate; if the combination satisfies m£>o 7 ( 7 ) — m D o < 150MeV/c 2 and 
cos 6 by 1 < 1-1, where Y' stands for the D°^{ r y)£ system, the photons are considered as a part of 
the tag system and removed from the signal system. 

At this point, we require that the event contains no charged tracks besides the B~ — > D°£~17(X) 
and B + — > -k®£ + v candidates (T extra = 0). We further require the total residual neutral energy of 
the event in the center-of-mass system to be less than 300 MeV (E ex t ra < 300 MeV). 
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To simplify the signal extraction procedure, only events with a sing le B+ -> ir°£ + i> candidate 
passing all the selections are kept. Only 1.6% of the MC signal events have more than one such 
candidate. 

Distributions of the selected events in cos B 7T oi and the difference between the reconstructed 
and fitted D° mass are shown in Figured These variables are used to extract the signal yield. 

3.2 Signal Yield 

We extract the total signal yield as 

N SB _ N SB 

N = N on -N oS -N MC atSB > ( 4 ) 

iV MC 

where N on , JV fj, and Nmc are the numbers of events in the signal region in the on-resonance 
data, off-resonance data, and simulated BB background, respectively, after the Z)°-mass sideband 
subtraction. Note that this subtraction is performed with the assumption that the D" mass non- 
peaking background is linearly distributed over the range ±10<r around the nominal D° mass. The 
good fit of a straight line to the sideband distribution indicates that this assumption is reasonable. 
The numbers N q q and Nmc are scaled according to the integrated luminosity of the on-resonance 
data. The "SB" in superscript indicates a cos# B7r o^ sideband region defined as —20 < cos9 BlT 0£ < 
— 1.5. The ratio (N^ — N^)/N M q is used to correct for the effect of data/MC discrepancies 
in the D°-mass peaking background yield selected in the signal region —1.1 < cos6 B7T ag < 1.0. 
We assume that this ratio has the same value in the signal region since both regions are mostly 
populated by events in which the tag-side decay is B~ — > D^°£~"P. The ratio is found to be 
1.02 ±0.11 ±0.35 where the first error is statistical and the second error is a systematic uncertainty 
derived from the variation of the ratio when it is computed in three statistically independent bins 
of the —20 < cos#B ]7r 0£ < —1.5 region. 

We find no off-resonance events in the signal and D° mass sideband regions. The statistical 
uncertainty on these numbers, derived from the corresponding Poisson intervals, is large. Moreover, 
this uncertainty would have to be scaled by a large factor (~ 8) to match the integrated luminosity 
of the on-resonance data. This would result in a large error on the background subtraction, and calls 
for a special treatment. Thus, we consider the alternative of computing the continuum background 
yield using a factorized efficiency (€mult) equal to the product of the individual efficiencies of 
the last two cuts of the event selection, T ex t ra = and E ex t ra < 0.300 GeV, which have a high 
power of background suppression. The direct use of emult would be justified if these two cuts 
were uncorrelated. The variation between the usual efficiency and €mult has been studied as the 
Textra and E ex t ra cuts were progressively relaxed. From the difference in the behavior of these two 
estimators, a correction factor has been derived and a systematic uncertainty has been added to 
the multiplied efficiency: Cmult*^ = (0.35±0.35)eMULT- The resulting off-resonance event yield in 
the signal region, after the Z)°-mass sideband subtraction and scaled to the integrated luminosity 
of the on-resonance data, is —0.3 ± 2.5. 

We find Nan = 52.0 ± 8.1 events in the signal region, —1.1 < cos9 BlT 0£ < +1.0, after the IP- 
mass sideband subtraction, of which 6.7 ± 2.0, 0.6 ± 0.9, and —0.3 ± 2.5 events are estimated to be 
B + B~, B°B° and non- BB backgrounds, respectively. The errors are statistical. The signal yield 
is therefore 44.9 ± 8.7 events. The corresponding yields before the L>°-mass sideband subtraction 
are 61 ± 7.8, 9.0 ± 1.9, 1.6 ± 0.8 and 1.5 ± 1.7 for on-resonance data, B+B~, B°B° and non-BB 
backgrounds, respectively. Figure |3] illustrates the yields for both data and MC samples and thus 
the significance of the signal. 
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Figure 2: Event distributions after all cuts (except on the variables shown), in the D mass 
- cos^y^ plane. The signal region is bounded by the black box: the reconstructed D° mass 
values must be within 3a of their fitted mean (a is the experimental resolution on the D° meson 
mass) and values of cos6 B lT oi must lie between -1.1 and +1.0. The —20 < cos^ Bj7r 0£ < —1.5 
region is defined as the cos 6g n oi sideband while the D°-mass sideband regions correspond to 
3cr < |m^o° — iri^o ed \ < 10<7. The MC BB and off-resonance distributions are not scaled to the 
integrated luminosity of the on-resonance data. 
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Figure 3: cosOs^e distribution, after all other cuts, for events having the reconstructed D° mass 
values within 3a of their fitted mean. The points are the on-resonance data and the error bars 
are statistical only. In order to illustrate the final yields and the significance of the signal, each 
distribution (on-resonance data, MC signal, MC B+B~ and MC B°B°) is rescaled to the number 
of events found after the D°-mass sideband subtraction. Note that the off-resonance data, due to 
their special treatment (see Section l3.2(l . are not shown. The final cos^g^c^ selection is between 
-1.1 and 1.0, delimited by the arrows. 
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3.3 Signal Efficiencies 

The signal yield is expressed in terms of the efficiency e data 7 computed using a sum of B + — > ir e + v 
and B + — > ir°n + v decays, as 

iV = 2e data 23(5+ -» 7r ^ + i/)7V B , (5) 

where the factor of two comes from using electrons or muons in the branching fraction, and Nb 
is the number of mesons in the data sample. Using the T(45) — » B°B° branching fraction 
Zoo = §> the number Nb equals 2(1 — foo)N BB -, where N BB - is the number of BB events in the data 
sample and the factor of two comes from having two B mesons in each event. 

We use the signal Monte Carlo simulation to provide an estimate of the total efficiency: 

e MC = (1.40 ±0.10) x 10 -3 , (6) 

The uncertainty is due to Monte Carlo statistics. 

We evaluate the data-MC difference of the B~ — > D°£~V(X) selection efficiencies using the 
double-tag events, in which both B mesons decay to D°£v(X). Properties of the B~ -> D°e~u(X) 
tags such as the composition of the D° decay channels are similar between the tagged-signal and 
double-tag events. The number of double-tag events is proportional to the square of the tagging 
efficiency after subtracting the small contribution from background. 

The selection criteria for the double-tag events follow the main analysis as closely as possible. 
In each event, we look for two D°£u(X) tags that do not overlap with each other. We remove all 
particles that are used in the tags and require that there be no charged tracks. The E ex t ra cut 
is increased from 300 MeV to 600 MeV; this reflects the fact that the residual neutral energy in a 
signal event comes from the tag-side B decay, and a double-tag event contains two such decays. 

From the study of the double-tag events we extract the efficiency correction factor 

pdata 

= 1.005 ± 0.039. 

The error includes both the statistical and systematic uncertainties. For the latter, we consid- 
ered the difference in the results when the selection criteria are varied, the residual background 
after the D°-mass sideband subtraction, due to possible non-linearities in the background vs. the 
reconstructed Dz mass, and the uncertainties in the exclusive B — > X c lv branching fractions. 

4 SYSTEMATIC UNCERTAINTIES 

The significant sources of systematic uncertainties considered and their impact on the measured 
branching fraction are summarized in Table ^ 

The systematic error is dominated by the size of the available MC samples. This includes both 
the signal MC sample for evaluating the signal efficiency and the generic BB MC samples for the 
background subtraction. 

The systematic uncertainties associated with the signal ir° efficiency, the identification of the 
signal lepton and its tracking efficiency are derived from detailed studies of the discrepancies in 
particle reconstruction between MC simulation and data. In particular, the systematic uncertainty 
in signal ir° efficiency is obtained by comparing the ratio of the t~ — » h~ yield to the 

r~ — > h~TT u yield in data and simulation, using e + e~ — ► t + t~ events. The uncertainty associated 
with the tagging efficiency was evaluated by using double-tag events (see Section f3.3|) . 
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Table 1: Fractional systematic errors on the branching fraction. 



Systematics 


a B /B 


MC BB statistics 


±4.6% 


ML signal statistics 


±7.3% 


(N* n B - NjgyN&f, ratio 


±6.0% 


Signal lepton tracking 


±0.8% 


Signal lepton ID 


±3.3% 


CI* 1 £32 

Signal 7r emciency 


±5.0% 


Tagging emciency 


±3.9% 


5 -» vr^ FF 


±1.5% 


B(B —>■ X c / U iv) background 


±1.1% 




±1.1% 


Total 


±12.9% 



For the B — > irlv form factor, we used the ISGW2 model 0]. The related systematic uncertainty 
is evaluated by reweighting ^3] the simulated signal events to reproduce the q 2 distribution pre- 
dicted by the Ball-Zwicky, HPQCD and FNAL calculations ^[^HHl- This uncertainty corresponds 
to the spread between calculations. The form factor affects the branching fraction through the q 2 
dependence of the signal emciency. Only the shape and not the normalization of the form factor is 
relevant for the measurement of the branching fraction. 

Background events are mostly from B — > X c lv decays, with a small contribution from B — > 
X u tv decays. For many of these decays, the branching fractions are not well known. To evaluate the 
associated uncertainty, the branching fractions of these decays were varied simultaneously within 
their uncertainties ^1]. The resulting variation in the B + — > ifil^v branching fraction was taken 
to be the corresponding systematic uncertainty. 

In addition to the studies discussed above, we performed crosschecks to ensure that the simu- 
lation was adequately modeling the data close to the signal region, before unblinding the data. In 
particular, we've studied the data-MC agreement in sidebands of cos 9 BlT 0£ and E extra , and observed 
no variation beyond expected statistical fluctuations. 

5 RESULTS AND SUMMARY 

Using event samples tagged by B~ — > D°£~V(X) decays, we determine the exclusive branching 
fraction B{B + — > tt°£ + v). From the signal yield and the efficiencies evaluated in Section |H1 we 
extract a preliminary result 

B(B+ 7r°£+u) = (1.80 ± 0.37 stat . ± 0.23 syst .) x 10" 4 . 

Table 12 summarizes the measurements of B(B — > niv) by the BABAR collaboration. Assuming 
isospin symmetry and the ratio of B lifetimes t^+/t^o = 1.081 ± 0.015 |14j . the measurements 
agree with each other with x 2 = 10.3 for 5 degrees of freedom, which corresponds to a one-sided 
probability of 7%. 

We plan to extract IV^I in the next version of this analysis, which will also include more data. 
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Table 2: Preliminary BABAR measurements of B{B — > irlv). The last row shows this measurement 
and the B° — > it^t^v measurement in Ref. [D] 

Technique B(B° -> n-j+v) x 10 4 B{B + -> 7r°£+^) x 10 4 

Neutrino reco. 1.41 ± 0.17 ± 0.20 0.70 ± 0.10 ± 0.10 

Hadronic tag 0.89 ± 0.34 ± 0.12 0.91 ± 0.28 ± 0.14 

Semileptonic tag 1.03 ± 0.25 ± 0.13 1.80 ± 0.37 ± 0.23 
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